RV survey of low-mass companions to sdB stars 
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Abstract. The origin of subdwarf B (sdB) stars is not fully understood yet since it requires high 
mass loss at the red giant stage. SdBs in close binary systems are formed via common envelope 
ejection, but the origin of apparently single sdB stars remains unclear. Substellar companions may 
be able to trigger common envelope ejection and help forming sdBs that appear to be single. 
Using a sample of high resolution spectra we aim at detecting small radial velocity (RV) shifts 
caused by such low mass (sub-)stellar companions. The RVs are measured with high accuracy using 
sharp metal lines. Our goal is to test the theoretical predictions and put constraints on the population 
of the lowest mass companions to sdB stars. 
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INTRODUCTION 

Subluminous stars are considered to be core helium-burning objects with a thin hy- 
drogen shell. In a Hertzsprung-Russell diagram these stars are situated on the Extreme 
Horizontal Branch (EHB). The evolution of these stars is far from understood since the 
formation of an sdB requires extreme mass loss at the tip of the RGB. The mechanism 
leading to this mass loss is still under debate [3, 4]. 

Formation scenarios can be roughly divided into single star and binary channels. 
About 50% of the known sdBs are found to be in close binaries with periods of few hours 
to several days consistent with the predictions of the common envelope (CE) ejection 
scenario [2]. Two main sequence stars evolve in a binary system. The more massive one 
becomes a red giant first and fills its Roche Lobe. Unstable mass transfer leads to the 
formation of a CE. The core of the red giant and the immersed secondary experience a 
loss of orbital energy which is deposited in the CE causing a spiral-in towards the center 
of mass. As soon as enough energy has been deposited in the envelope, it is ejected. 
The inwards migration stops and the stars end up in a close binary. This scenario is well 
suited to explain the known population of close sdB binaries. 

But how are the apparently single sdBs formed? Soker [10] suggested that substellar 
companions like planets or brown dwarfs may be able to trigger CE ejection as well. 
In order to survive the common envelope and avoid evaporation inside the envelope the 
secondary should have a minimum mass of more than lOMj. 

A similar scenario has been proposed for the formation undermassive single white 
dwarfs by Nelemans and Tauris [7]. 

Several substellar objects have recently been detected around hot subdwarfs in wide 
orbits, showing that planets can survive the red giant phase of their host star [9, 5,8]. 
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FIGURE 1. The plot displays a section of a FEROS spectram of HD 205805 featuring exemplary sharp 
metal absorption Unes. 



Recently, Geier et al. [1] discovered a candidate substellar companion in close orbit 
around the bright sdB HD 149382, which perfectly fits into the pattern predicted by 
Soker. Systems like that may be quite common and could have remained unnoticed up 
to now. 

The goal of our project is to put constraints on the population of the lowest mass 
companions to sdB stars by measuring accurate radial velocities from multi-epoch high 
resolution spectra. 



OBSERVATIONS AND DATA REDUCTION 

Up to now our survey consists of 23 bright sdB stars with visual magnitudes ranging 
from 10 to 14. We used high resolution spectra obtained with the FEROS instrument 
mounted at the ESO/MPG-2.2m telescope. The spectra have a resolution of 48000. 
Spectra were reduced, calibrated and corrected for earth's orbital motion with the MI- 
DAS package using the FEROS reduction pipeline. 

All programme stars are single-lined objects without visible spectral features of a 
companion. Furthermore, those binaries with high radial velocity variability have been 
excluded because their unseen companions must be quite massive. Each star has been 
observed several times (from 2 to 15) on timescales ranging from one day to some years. 
As an example we display a section of the FEROS spectrum of HD 205805 in Fig 1. 
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FIGURE 2. Spectral features of HD 205805 fitted with FITSB2. The left hand panel displays metal 
absorption lines suitable for RV measurements, in the right hand panel exemplary night-sky emission 
lines used to check the wavelength caUbration are shown. The latter are shifted in RV due to heliocentric 
correction applied to the spectra. 



MEASUREMENT 

Subdwarf B stars are hot (Teff — 20000 — 35 000 K) and their atmospheres radiative. The 
dominant features in their spectra are the hydrogen B aimer and several helium lines. 
But the spectra also show up to ~ 50 weak metal absorption lines which are very sharp. 
These features are best suited to measure radial velocities. 

The atmospheres of sdB stars are considered to be stable. The only known effect 
other than a close companion that may lead to RV variability are pulsations. Pulsating 
sdB stars do indeed exist with periods between 90s and 3h and RV amplitudes of few 
km s^^ at most [4]. 

For our analysis we chose a set of sharp, unblended metal lines situated between 3600 
and 6600 A (see Fig. 2, left hand panel). Accurate rest wavelengths were taken from the 
NIST database. Depending on the atmospheric parameters of the star and the quality of 
the data a subset was used. 

The RV measurement consisted of two steps. First, we fitted Gaussian and Lorentzian 
profiles to each line separately using the SR\S routine (Hirsch priv. comm.). The re- 
sulting mean RVs had standard deviations ranging from 0.8 to 2.7kms~^ Second, we 



used the FITSB2 routine [6] to perfomi a simultaneous fit of all suitable lines. Uncer- 
tainties were calculated using a bootstrapping algorithm and found to range from 0.1 to 
0.5kms~\ while both methods result in consistent mean RVs. 

The scatter of individual RV values measured with SRA.S around the mean value was 
reviewed for possible systematic effects like correlations with certain elements. No such 
effects could be found. To check the wavelength calibration for systematic errors we 
used telluric features as well as night-sky emission lines (see Fig. 2, right hand panel). 
Having their origin on earth these lines should have zero RV. The FEROS instrument 
turned out to be very stable. Usually corrections of less than 0.5kms^^ had to be applied. 

In order to derive the degree of precision we can achieve with our method, we gen- 
erated synthetic model spectra with realistic atmospheric parameters including Poisson 
noise with S/N — 1000 and determined the RV in the way described above. The un- 
certainties turned out to range from 0.1 to 0.4kms^^. The accuracy is limited by the 
resolution of the spectrograph, the number of features and the intrinsic thermal broad- 
ening of the lines. 

RESULTS 
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FIGURE 3. The diagram shows the result of RV determination for two spectra of EC 20106—5248 taken 
6 months apart. Each measuring point represents a metal absorption line. The horizontal lines marking 
mean RV values of the spectra indicate a significant shift in radial velocity. 



TABLE 1. Summary of RV results for our sample obtained using two different algorithms. 
RVs for each line were measured independently and averaged with SR^S (left hand colunm) 
whereas FITSB2 uses a simultaneous fit (right hand column). 
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Four stars were found to be significantly (another two marginally) RV variable. The 
shift in RV ranges from 0.5 to 12.7kms^^ (see Tab. 1). As an example we display 
two RV measurements for EC 20106—5248 in Fig. 3. Follow-up photometry and high 
resolution spectroscopy are necessary to exclude pulsational variability and derive the 
orbital parameters of these binaries. Constraints can then be put on the companion 
masses. Even if there are no low mass companions we would expect to find a certain 
fraction of systems with small RV shifts caused by more massive secondaries seen at 
low inclination. A larger sample is needed to exclude this case. 

No RV variations were found in the rest of the sample. Nevertheless, the errors allow 
to derive upper limits for Mcomp sin / of undetected substellar companions in close orbits 
from the binary mass function (Eq. 1). 

M^oj^p sin^ i PE? 

^^^^ ^ (Mcomp+M,dB)2 " 2^ 

Soker [10] suggested that these objects should be more massive than lOMj and should 
have orbital periods shorter than 10 d. Adopting this maximum period and assuming a 
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FIGURE 4. The histogram shows upper mass Umits for unseen companions in close orbits around sdBs, 
where no significant RV shifts could be detected (left hand axis). These are derived from upper Umits for 
the RV semi- amplitudes (right hand axis). 



canonical mass of 0.47 Mq for the sdB, the errors of the RV measurements can be used 
as upper limits for the RV semi-amplitudes (K). 

As can be seen in Fig. 4, tight constraints can be put on possible substellar compan- 
ions, if the measurement error is smaller than ~ 1 .0 km s~ ^ . Given data of perfect quality, 
close-in planets with masses as low as 0.3Mj can be firmly excluded if our assumption 
about the orbital period is correct. 

In conclusion ~ 26% of the apparently single stars in our sample show small RV 
variations. The true fraction is expected to be higher, because the accuracy is limited 
by the S/N in most cases. Companions with masses higher than lOMj in close orbits 
(P < lOd) can be excluded for 61% of our sample. 

One possible explanation for the fractions found is that the assumed CE scenario is 
not the only channel for sdB formation. However, even with CE ejection being the only 
channel, the companion may be destroyed during the CE phase either by evaporation or 
a merger with the stellar core leaving behind a single sdB. Better statistics and reliable 
predictions from theoretical calculations are required to solve this problem. 
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